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ABSTRACT 

Aims. We describe a 6-12 GHz analogue correlator that has been developed for use in radio interferometers. 

Methods. We use a lag-correlator technique to synthesis eight complex spectral channels. Two schemes were considered for sampling the 
■ cross-correlation function, using either real or complex correlations, and we developed prototypes for both of them. We opted for the "add and 
square" detection scheme using Schottky diodes over the more commonly used active multipliers because the stability of the device is less 
critical. 

Results. We encountered an unexpected problem, in that there were errors in the lag spacings of up to ten percent of the unit spacing. To 
' overcome this, we developed a calibration method using astronomical sources which corrects the effects of the non-uniform sampling as well 
. as gain error and dispersion in the correlator. 
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1. Introduction 

High brightness sensitivity is essential for many observations 
in radio astronomy, for example, observations of the cosmic 
microwave background (CMB). Once the angular resolution 
required has been fixed by the size of the instrument in wave- 
lengths, the brightness sensitivity is determined by two factors: 
the system temperature and the observing bandwidth. System 
temperatures have been reduced dramatically by the use of 
cryogenic HEMT amplifiers, to the point where the atmosphere 
rather than the amplifier may be the dominant noise source in 
the system, and thus little further improvement is possible for 
ground-based instruments. Increased continuum bandwidth is 
therefore an important route to increased sensitivity. 

Radio interferometers with any significant bandwidth re- 
quire the band to be broken up in to sub-bands, both because 
spectral resolution may be intrinsically required for the obser- 
vation, and to overcome the effects of chromatic aberration. 
Digital correlators are commonly used since they provide the 
most convenient way of generating high spectral resolution and 
a large range of time delays. However, for very high bright- 
ness sensitivity continuum observations which demand max- 
imum possible bandwidth with only modest spectral resolu- 
tion, analogue correlators are competitive. Analogue correla- 



tors have been used in several recent interferometers dedicated 
to CMB observations dp' Sullivan et all 1 1995b ILo et al l 12001 



Padin et alJEool iLeitch et al.ll2002l: I Watson et alJl2003h . 

To achieve the low spectral resolution typically required 
(of order ten spectral channels, rather than the hundreds 
or thousands typical of digital correlators) there are two 
basic solutions. One is to split the broadband signal di- 
rectly into small er bands with filte r banks or multiple down- 
converters (eg. iPadin et al . 2001). The other is to use a 
single broad intermediate frequency (IF) band and imple- 
ment a Fourier transform spectrometer, or lag correlator. 
This forms the cross-correlation between antennas at multi- 
ple delays, and the outputs are Fourier transformed to es- 
timate the cross-frequency spectrum. In this paper, we de- 
scribe the development of an analogue 6 GHz bandwidth 
lag correlator with eight complex bands. The correlator is 
currently operational in t he Arcminute Micro k elvin I mager 
(AMI; see, for example | Kneissl et al.1 l200ll: Uonesl 120021: 
Kaneko & the AMI Collaborationl 20061) . a new radio interf er- 
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ometer operating over the 12-18 GHz band. AMI was specif- 
ically designed to carry out a survey for clusters of galaxies 
through the Sunyaev-Zel'dovich effect ( Sunyaev & Zel dovichi 
Il972h . Broadband anal ogue correlators with s imilar concepts 
have been described b y lHarris & Zmuidzinasl d200ll) : iLi et al.l 
( 2004 ): [Roberts et aD (submitted). A feature that distinguishes 
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our design from these correlators is the detectors. Rather than 
active multipliers, we chose the "add and square" detection 
scheme using zero-biased Schottky diodes. They are simpler 
to operate because no bias current is required and their stabil- 
ity is not as critical as it is for active devices. It is technically 
challenging to design IF components that span over more than 
one octave in frequency so we chose an IF band of 6-12 GHz. 
To reduce costs and complexity, we implemented the correlator 
using microstrip technology and each baseline was fabricated 
on a single substrate. 

Section [2] of this paper discusses the basic relationships for 
a broadband lag correlator and discusses real vs complex corre- 
lators and phase- switched multipliers. Section [3] describes the 
detector circuit in our design. Section [4] then describes the over- 
all hardware implementation, and Sec.[5]the calibration scheme 
implemented to recover accurate spectra from the lag data in 
the presence of dispersion and lag spacing errors. Section [6] 
describes issues related to the manufacturing and practical im- 
plementation of the correlator, with conclusions in Sec. [71 

2. Correlator Overview 

2.1. Fourier Transform Correlator 

The response of an interferometer baseline with a single real 
correlator having a rectangular passband of width Ay is given 
by 



sin[7rAv(r g - Ti)] 

7TAv(Tg - Ti) 



■ COS 



27r(v L oTg - Vi F (T g - Ti)) + L( 



(1) 



where Ro is a normalising factor, T g and Ti are geometric and 
instrument delays (i.e. the delays introduced at the radio fre- 
quency (RF) and IF respectively), v LO is the local oscillator 
(LO) frequency, v IF is the IF frequency and LO is the phase 
difference between the LO at each aerial. The cosine term de- 
scribes the fringes that change with changing delay. The path 
compensator inserts an instrument delay Ti which (approxi- 
mately) balances out the geometrical delay. The sine term pro- 
duces a fringe envelope, which maximises the response if the 
residual delay, (r g - Ti), is small. For simplicity, the phase dif- 
ference in the LO at the two mixers, LO , can be dropped as it 
is easily removed by calibration. Substituting for the LO fre- 
quency v LO using the relationship v LO - Vi F = v RF for lower side- 
band reception, 



sin[7rAv(T g - Ti)] 
R = R — — cos 

7TAv(Tg - Ti) 



2tt( 



v rf t 2 + v if t 



(2) 



Chromatic aberration is a problem because of the extent of 
the observing field of view. While it is possible to introduce the 
correct instrumental delay to compensate for the geometrical 
delay at the centre of the field, a radio source a small angle AO 
away from the field centre will have an additional geometric 
delay DA6/c, where D is the baseline length projected on to 
the plane perpendicular to the source direction. If it is required 
to image with minimal loss of sensitivity to the edge of the 
primary beam of an interferometer with an antenna diameter d 

Av/vr 



(see eg. iThompson et al. 1200 lh . To satisfy this requirement 
for an instrument with centre frequency v RF = 15 GHz and a 
maximum ratio of baseline to dish diameter D/d ^ 10 we 
chose to split the full 6 GHz bandwidth into eight sub-bands 
of 0.75 GHz each. 

The basic layout of the Lag or Fourier transform correlator 
is illustrated in Fig.Hk. The signals A(t) and B(t) from the two 
antennas of a baseline are correlated at discrete time delays Ti 
with step size Sr. The set of correlated signals form the cross- 
correlation function: 



R( Tl ) 



A(t) B(t + T { )&t, 



(4) 



where T is the integration time. This type of correlator mea- 
sures the temporal coherence function of a signal, where the 
coherence function is an indirect measure of the signal's fre- 
quency spectrum. The cross-power spectrum (hereafter abbre- 
viated to the spectrum) can be recovered by applying the dis- 
crete Fourier transform (DFT) to the cross-correlation function; 



Ti 



■2njv k Ti 



(5) 



d/D 



« 1 
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Figure [2] illustrates typical cross-correlation functions of a 
flat- spectrum point source at the telescope's phase centre for a 
narrowband and broadband system with rectangular passband. 
The function is given by Eq.[2]with T g = and Ti as the vari- 
able. The cosine term has a period determined by the centre IF 
frequency v IF and the sine envelope is defined by the bandwidth 
Ay. As the source moves across the sky, the cosine pattern drifts 
in instrument delay space. The envelope encodes the spectral 
information and its Fourier transform gives the spectrum of 
the observed signal. The desired information is contained in 
the envelope and it must be sampled at Nyquist sampling step, 
6t n = l/(2Av). 

A point source at the edge of the field of view will introduce 
an additional geometrical time delay AT g into one arm of the 
correlator. The additional delay will offset the cross-correlation 
function away from the centre lag. We define a narrowband 
correlator to be one where this delay AT g is shorter than the 
Nyquist sampling step Sr^. In order to avoid chromatic aberra- 
tion, no additional lags are necessary. However in a broadband 
system, where AT g is greater than Stn, the offset cannot be ne- 
glected. For our design parameters, AT g can be up to ±6 Stn 
so the cross-correlation function has to be sampled over a wide 
interval. We sample 16 lags at -8£tn, -7Stn, • • • +7Stn- After 
applying the DFT, this splits the passband into eight complex 
sub-bands of 0.75 GHz each. Most of the signal power is con- 
centrated in the centre of the cross-correlation function's en- 
velope. Its position in instrument delay space depends on the 
position of the source in the field of view and must be within 
the range of the sampled lags. Figure[3]shows simulations of the 
amplitude response for each output channel of a 16-lag correla- 
tor as a function of the centre position of the cross-correlation 
function. By using 16 lags for the correlator, its sensitivity is 
not significantly reduced for sources within the telescope's field 
of view. 
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Fig. 2. The cross-correlation functions for a broadband (v IF = 
9 GHz, Ay = 6 GHz) and a narrowband system (Ay = 
0.75 GHz). The horizontal scales are instrument delay in elec- 
trical length. The gray envelopes are the bandwidth patterns 
and the open circles are the lag positions. A source at the edge 
of the field of view will introduce an additional geometrical 
time delay Ar g and offset the cross-correlation function over 
the range of the arrow. 



Fig. 1. (a) The lag correlator cross-correlates the signals from 
a pair of antennas at discrete delay steps (denoted by St). 
Applying a discrete Fourier transform to the set of lag data 
gives the complex power spectrum. The real correlator samples 
at 16 lags, (b) The cross-correlation function and the spectrum. 
The gray envelope is the bandwidth pattern and the open cir- 
cles are the positions of the lags. For a source at the centre of 
the field of view and no geometric delay, the imaginary com- 
ponent of the spectrum is zero. In general, the real correlator 
has complex spectra at positive and negative frequencies. 

2.2. Signal recovery 

A narrowband signal of centre frequency v can be sampled by 
a simple single-lag correlator. The signal received at each an- 
tenna is a real quantity but it is convenient to express it in com- 
plex notation. The signals at the two antennas are the real com- 
ponents of 

A(t) = a e K2nvt+ ^ and B(t) = b e j2nvt , 

where is the signal phase which depends on the position of 
the source and geometry of the interferometer. After correlating 
these signals, the output of the correlator is 

Re[A(0 B*(t)] = Re[ab e j(f) ] = ab cos 0. (6) 

We take the complex conjugate of B(t) because we are inter- 
ested in the product of the real parts. The quantities we want to 
measure are the amplitude ab and the phase 0. To measure them 
both simultaneously, we would need two measurements. We 
could make the second measurement by inserting a 90° phase 
shift in one arm so that 

A'(t) = ae j(2nvt+(p+7T/2) . (7) 



1.4 




Position of peak from lag centre /lag 

Fig. 3. The amplitude response of each frequency channel 
of the correlator to the broadband cross-correlation functions 
(Fig. [2) positioned away from the lag centre. This is equivalent 
to point sources away from the pointing centre. The sensitivity 
of the correlator is fairly constant over the field of view, cor- 
responding to a maximum of ±6 lags. The structures and the 
reduced sensitivity in the edge sub-bands 1 & 8 will be dis- 
cussed in Sec. 15.21 

The second measurement will be 

Re[A r (0 B*(t)] = Re[ab e M+nl2) ] = ab sin 0. (8) 

From these in-phase (equation 3) and quadrature (equation 4) 
components, we can simultaneously measure the amplitude and 
phase of the signal. As mentioned earlier, the Fourier trans- 
formed output of the lag correlator is a complex spectrum. 
This is a way of expressing the amplitude and phase of each 
sub-band. Similarly, we could express the amplitude and phase 
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measured by the single-lag correlator as a complex value. We 
stress that the cross-correlation function is a real quantity and 
the complex notation is for mathematical convenience. To re- 
flect this notation, we will refer to the in-phase and quadrature 
components as real and imaginary. 

In a multiple-lag correlator, however, we have two options. 
We could either sample just the in-phase components of the 
cross-correlation function (real correlator - see Fig.[Tk) or both 
the in-phase and quadrature components (complex correlator - 
see Fig. Hk). The choice will affect the form of the spectrum 
S (vk). But the information we recover is equivalent. 

2.2.1. Real correlator 

Here, In real components of the cross-correlation function are 
sampled at Nyquist rate with step size Stn ; an electrical length 
of 25 mm or an equivalent sampling frequency of 12 GHz for 
a 6-GHz bandwidth signal. The cross-correlation function is 
sampled from -n • Stm to +(n - 1 ) • Stn, where n is the number of 
complex sub-bands [J. The symmetry between the positive and 
negative halves of the spectrum means that half the information 
in the spectrum is redundant. 16 independent measurements of 
the cross-correlation function give 8 independent complex sub- 
bands (see Fig.[U)). 

2.2.2. Complex correlator 

The second method samples both the real and imaginary com- 
ponents of the cross correlation function but at only half- 
Nyquist rate (6 GHz sampling frequency, or an equivalent elec- 
trical length of 50 mm). We express the measurements as a set 
of n complex numbers. The step size is 2Stn and samples span 
from -n • Stn to +(n - 2) • Stn- The resulting spectrum contains 
only positive frequencies (Fig.0]3). To measure the imaginary 
component, one arm of input signal has a broadband 90° phase 
shiftei@. 8 complex measurements will give 8 complex sub- 
bands. The single- sided spectrum explains why half-Nyquist 
sampling is sufficient. It can also be understood in terms of 
information content: The information in 8 complex measure- 
ments is the same is 16 real measurements. 

To explore possible practical differences between the two 
methods, we built prototypes for both schemes. 

2.3. Cross-Correlation and Phase-switching 

There are several ways of correlating two analogue signals. In 
a single baseline interferometer, two antennas measure random 
noise-like signal voltages A(t) and B(t). We will represent these 
by random variables a and b. Active multipliers will directly 
give the product ab. Active components have higher sensitivity 



1 Generally, the spectrum will be complex in both positive and 
negative frequencies. These two components are Hermitian; S(v) = 
S*(-v), where the asterisk denotes the complex conjugate. The rela- 
tionship could be observed from the symmetry in the spectrum. 

2 The phase shift produces a Hilbert transform of the original sig- 
nal; for a single-sided spectrum, the real part of the cross-correlation 
function is the Hilbert transform of its imaginary part. 
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Fig. 4. (a) The complex correlator samples at only eight lags but 
with two measurements per lag and twice the delay steps. Using 
the 90° phase shifters, it samples the in-phase and quadrature 
components of the cross-correlation function at each lag. (b) 
The cross-correlation function and the spectrum. In general, the 
spectrum recovered by the complex correlator is single-sided. 

but their stability is often more c ritical. The re lative merits of 
active multipliers are discussed in lHarrisI ( 2003 ). 



We use a very simple concept, described by lRyld (I1952|) as 
"add and square" detection. The two signals from the antennas 
are summed to give a + b. The summed signal is then passed 
through a non-linear device, in our case a diode, with an output 
of the form 



(a + b) 2 = a 2 + 2ab + b 2 



(9) 



The total power terms of each antenna a 2 and b 2 dominate the 
output, and the cross-correlation term ab needs to be extracted 
before recording the data. We use phase- switching to remove 
the total power terms, which also reduces slowly varying offsets 
and cross-talk in the system. The detection scheme is illustrated 
in Fig. [3 The two signals are modulated by Walsh functions / 
and g. These are periodic bi- valued functions which have the 
property that all the functions in a given family are orthogonal 
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to each other, and to their product with any other Walsh func- 
tion. Phase shifting the signal by 180° applies the modulation. 
The output from the correlator will be 

(/a + gb) 2 = /V + 2/gab + g 2 b 2 . (10) 

The products of Walsh functions with themselves, f 2 and g 2 , 
are constants. The product of two different Walsh functions, fg, 
is another Walsh function. To extract the cross-term, the output 
is demodulated by fg and integrated over the repeat period T w 
of the Walsh function: 

1 r^ w 1 r^ w 

fg A 2 + 2/gab + g 2 b 2 ) dt oc — ab At. (1 1) 

Jo Jo 

Phase- switching also removes any slowly varying drift on time- 
scales longer than the repeat period. It also suppresses cross- 
talk because signals from other antennas will be modulated 
with an orthogonal Walsh function. 

To achieve the full signal-to-noise ratio (SNR) in such a 
correlator, it is also necessary to form the "subtract and square" 
product (a - b) 2 . This can be seen by considering that each 
voltage being correlated consists of a signal term which is cor- 
related between the two antennas and noise term that is uncor- 
rected. Since the noise terms are uncorrelated, two orthogonal 
signal combinations can be made from them with uncorrelated 
noises. The noise on the sum signal, a + b, is thus uncorrelated 
with the noise on the difference signal, a - b. Each can then be 
squared to form estimates of the product ab with independent 
noises, and these can be added to form the full signal-to-noise 
product. A full plus-minus system is thus equivalent to a sin- 
gle real multiplication, and two such multipliers would make a 
single complex correlator. 

A second correlator is thus equipped with a 1 80° phase shift 
in one arm. The output from this second correlator will be 



(/a-gb) 2 



fa 



■2/gab + g 2 b. 



(12) 



In principle the plus and minus correlator outputs can be dif- 
ferenced in hardware and the signal then demodulated and read 
out. This cuts down the necessary readout electronics by half. 
However, our final system demodulates and reads out the out- 
puts from each correlator separately. This enables individual 
errors in the correlator lag positions to be corrected before the 
signals are combined, (see Sec. 14.21) . 

The full schematic drawings of the real and complex corre- 
lators are shown in Figs. [6^ and[6j). The components for these 
correlators include broadband (6-12 GHz) signal splitters, slot- 
line 0° and 180° phase shifters, 90° microstrip phase shifters for 
the complex correlator and the detector circuits, including the 
diodes. Our desig ns of the splitte r and the phase shifters have 
been described in lHoller & Jonesl (submitted). 

3. Detector 

The detector circuit is the core of the correlator. We use a 
diode power detector circuit (Fig. [7]) comprising a Schottky 
diode operated in the square law region. For V in «c W, where 
Vt = k B To/e ^ 26 mV is the thermal voltage at room tem- 
perature 7p, the output vo ltage is proportional to input power 
faewlett-Packar dllAN 986h . 






Inte- 
grator 


) — - 







<ab> 



8(/ 2 a 2 +g 2 tf)+/^b 



Fig. 5. The "add and square" analogue correlator including 
phase switching. The signal from the detector is amplified and 
most of the total power terms are removed by the highpass fil- 
ter (HPF). The HPF attenuates the direct current (DC) signal 
by a factor e. Any remaining total power and DC offsets are 
removed by phase- switch demodulation. (•) is the expectation. 



For very small input signals a bias current could be used 
to increase the output. However, the output would be very sen- 
sitive to changes in this current with, for example, tempera- 
ture, so we chose to use zero-biased diodes. The disadvantage 
is that they require a relatively high input power in the region 
of -13 dBm or 0.05 mW. But they are very stable, as long as 
they are kept at a constant temperature. 

Tradeoffs had to be made between the most important per- 
formance measures of the detector. These include sensitivity 
over the full band, stability, the video bandwidth (the lowpass- 
filtered bandwidth after the detector), avoiding it behaving like 
a peak detector and frequency matching. 

The sensitivity largely depends on the value of the load re- 
sistor R L . A high value gives a higher output voltage for the 
same input power, but also reduces the video bandwidth for the 
signal behind the diode. Phase- switch modulation introduces 
sharp transitions in the signal. Slews over the transition will de- 
grade the orthogonality between the signals. For phase- switch 
signals with a maximum frequency of 1 kHz, we required a 
minimum video bandwidth of about 1 MHz. The bypass capac- 
itor also has a large influence on the video bandwidth, but its 
value should not be too low since it acts as a ground for the RF. 

Under certain circumstances, the detector could operate as 
a peak detector, where a large voltage peak charges up the by- 
pass capacitor Cb and no more current can flow through the 
diode. This can be avoided if the detector is operated in the 
small signal regime where the output voltage is comparable to 
the thermal voltage Vt. In this regime, the detector circuit will 
act as an integrator. Additionally, when operated in that regime, 
the detector's output voltage is proportional to the input power. 

After testing several circuits the diode MSS20-146 from 
Metelics was chosen together with R L = 10 kQ and Cb = 10 pF. 
The results can be seen in Fig. [U The passband is relatively 
flat and rather upwards sloping, which is advantageous, since 
most other microwave components will add a slope in the op- 
posite direction. The features in the passband are very simi- 
lar for all tested types of diodes, which means they can be at- 
tributed to the matching circuit rather than the diode. A plot 
of the diode response versus input power shows a square law 
region up to about 30 mV ^ Vt and is consistent with the sim- 
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(a) Real correlator schematic layout (b) Complex correlator schematic layout 



Fig. 6. The full schematic of the real (a) and the complex (b) correlators for a single baseline. The real correlator samples at 16 
lags in steps of 5tn, while the complex correlator samples at 8 lags in steps of 2Stn. At each lag of the complex correlator, the 
signal is split and one arm is phase-shifted by 90°. In both cases, there is a 180° board with one arm phase-shifted by 180°. 0° 
phase shifters of similar slotline design are used in corresponding positions to match the response with that of the 1 80° phase 
shifters. 



ulated model. The measured video bandwidth was in the range 
of 800-1 100 kHz. 




Fig. 7. The power detector circuit. The forward (direct) current 
through the diode produces a voltage across the load resistor 
R L . The bypass capacitor Cb looks like a short for the RF signal. 



so the impedance is complex as well as depending on input 
frequency. The minimum absolute value of the impedance for 
our frequency band can be calculated to about 100 £1 To keep 
reflections and signal scatterings between channels to a min- 
imum, we required S n to be low across the whole band. We 
used two element matching with a series microstrip line and 
a shunt line. In addition, the high impedance of the diode is 
reduced through a 100 D parallel lumped resistor at the cost 
of losing sensitivity. The resistor will bypass part of the signal 
away from the diode which will not be detected. The character- 
istics of the final detector circuit are shown in Fig. [9] 
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Fig. 8. (a) The detector response against power for a 9 GHz sig- 
nal. The square law region extends to about - 1 1 dBm input 
power, (b) The passband of the detector (including the match- 
ing circuit) from 5 to 13 GHz. In both plots, the solid lines are 
the measurements and the dashed lines are simulations. 



Matching the diode to the 50 Q microstrip circuit was an 
important task, both to maximise power transfer to the detec- 
tor and to avoid setting up standing waves in the correlator sig- 
nal distribution circuit. The diode includes reactive components 



S11 




Frequency (GHz) 



Fig. 9. A photograph of the final detector circuit (left) and the 
result of the matching (right). The plots show the measured 
input reflection (black curve) together with a simulation (gray 
curve) . 
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4. Correlator Hardware 



4.2. RF measurements 



4.1. Design 

In the final design, 0° and 180° correlations for a single base- 
line (see Fig. [6]) were implemented on two separate boards. 
Figures ITQk andfTQb show the physical layout of the 0° real and 
180° complex boards. The designs show all the components; 
splitters, phase shifters and the detector circuits. To combine 
the signals before the detectors, we used Wilkinson dividers in 
reverse. Each lag has different delay lengths. The layouts were 
designed to avoid right angle bends and to keep components at 
a distance to try to minimise possible coupling. The layout of 
the complex correlator was very similar to the real correlator. 
The most significant differences are the 90° phase shifters and 
different delay lengths. 



Phase shift 




Phase shift Output to Network Analyzer 

(a) Prototype 0° Real Correlator 



Phase shift 




Phase shift 



Output to Network Analyzer 



(b) Prototype 180° Complex Correlator 

Fig. 10. The physical layout of the 0° real correlator board (a) 
and the 180° complex correlator board (b). The corresponding 
180° real and 0° complex correlator boards are not shown. For 
these prototype versions, one detector was swapped for a 50 £1 
output to enable measurements on the network analyser. The 
numbers on the real correlator correspond to the lag orders. 



The transmission characteristics of the correlator through to the 
input to the diode detector have been plotted in Fig. [TU The 
signal has cascaded through the 1 80° slotline phase shifter, five 
signal splitters and one signal combiner (and one arm of the 
90° phase shifter in the case of the complex board). A lossless 
system would show a flat response at about -15 dBm. But mi- 
crostrip introduces loss, which is proportional to frequency and 
causes the slope across the band. The slope can be corrected 
for by inserting equalisers before the correlator input and is 
therefore of small concern to us. Large ripples across the band 
would be a problem, since they reduce the signal-to-noise per- 
formance of the correlator. 
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Fig. 11. Left: signal transmission through the correlator lags. 
Right: correlator input reflection. The results for the Complex 
and Real boards are very similar. 



The performance of the phase shifters inside the correla- 
tor is shown in Fig. [12] The 180° slotline phase shifter shows 
an exceptional behaviour with a maximum error of ±6° over 
the band of 5-13 GHz. The 90° microstrip phase shifter has a 
maximum error of about ±15°. 
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Fig. 12. (a) Performance of the 1 80° phase shifter in the system. 
For the performance of the 90° phase shifter (b), the difference 
between these two curves has to be taken. 



The frequency response, including all microstrip compo- 
nents, the detector circuit and the video band buffers of the pro- 
totype system is shown in Fig.[13la). There are some variations 
across the band but it is relatively flat with a gentle slope. The 
rising passband is similar to the detector's passband (Fig. [8J)) 
but there is a noticeable dip at 6 GHz. The structures in the 
passbands of the other lags differ slightly and the average sig- 
nal levels vary by a factor of 2 to 3. Therefore calibration of the 
correlator will be necessary. 
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Fig. 13. (a) Typical passband of a full lag including all RF com- 
ponents plus the detector circuit and the two first stages of am- 
plification. This shows a dip at 6 GHz compared to the response 
of the detector in Fig.[8j). (b) Dispersion in the outer most lags, 
which have the largest differential path length. The slope across 
the band is ±25°. 



Besides being lossy, particularly at higher frequencies, mi- 
crostrip also exhibits dispersion. For our substrate with e r = 9.7 
a physical length of one wavelength (1.02) at 6 GHz is 2.065/1 
at 12 GHz. However, this will have less influence on the cor- 
relator than one might expect. The two opposing signals in the 
correlator travel through similar length of micro strip line so 
dispersion is compensated. It is only the differential length of 
stripline that produces dispersion and the worst affected are the 
two outer lags. Figure \\3\ b) shows a measurement for one of 
these lags. It shows a phase slope of ±25° across the band of 
6-12 GHz, as expected from the distance to the zero lag. This 
implies that the electrical position of the point of correlation is 
a function of frequency. This effect will be compensated by the 
calibration process that will be discussed in Sec. 15. II 

The performances of the real and complex correlators were 
comparable. The ability to directly measure the amplitude and 
phase at each lag in the complex correlator is attractive for some 
commissioning tests. But calibrating for the small errors in the 
90° phase shifters over the passband adds an extra layer of com- 
plexity. From a practical view point, the more compact format 
of the real correlator fitted better into standard 6U racks. For 
these reasons, we chose the real correlators for the final sys- 
tem. 

One set of measurements did not turn out as expected. The 
physical distance between the lags varied much more than an- 
ticipated. As a consequence, the cross-correlation function is 
not uniformly sampled and a standard DFT cannot be used to 
recover the true spectrum. Instead, a means of calibrating the 
correlator had to be found. The exact positions of the correla- 
tions could be measured and this information could be used to 
estimate the true spectrum. The lag position errors on the pro- 
totype board were apparently random with values of 5-10% 
of the ideal lag spacing. However, in the final production se- 
ries, the lag errors were found to be smaller than in the pro- 
totype boards. Nevertheless, they will cause some reduction 
in the SNR. Such lag errors have also been found in other 



analogue bro adband correlators (IHarris & Zmuidzinasl 12001 
Roberts et al.1 submitted). 

The causes of the lag errors are unclear. The lag errors are 
not consistent between boards so a design error is ruled out. 
By replacing lumped components, we determined that the lag 
spacings are insensitive to the exact positions of the compo- 
nents. Inspection of the boards showed no problems with the 
printing or production of the board. A solder mask layer is 
printed to protect the lines and restrict solder flowing along the 
tracks. A variation in this layer could influence the effective di- 
electric constant. We manufactured a correlator board without 
a solder mask but the lag errors persisted. Before the thin gold 
layer can be plated on the tracks, the surface of the copper is 
roughened. It is conceivable that the roughened surface could 
increase the path length at high frequencies because most of 
the current flows along the surfaces of the lines. To test this, 
we also manufactured a correlator board with no gold coating. 
Reliable measurements for this board could not be made be- 
cause of poor contact to the connector. 

We found that groups of lags had similar lag errors to each 
other. This indicated that these groups may have been affected 
by a long arm of the delay line shared by the whole group. 
Variations in the relative permittivity of the dielectric e r could 
be a cause but the errors in neighbouring lags are too large to be 
explained by this. Variations in e r of ±0.5 would be necessary 
to explain the average errors and only if the dielectric varia- 
tions exactly followed the geometry of the board. This level 
of variation in e r is much higher than the manufacturing toler- 
ance. Whatever the cause may be, the lag spacings are stable 
over time. This enables us to correct for the lag errors using a 
calibration scheme we have developed for the broadband cor- 
relators. 

5. Calibration 

5.1. Method 

The eight complex frequency channels can be recovered from 
the 16 lag data by taking a direct Fourier transform (DFT). The 
DFT kernel can be conveniently expressed with the 16 by 16 
DFT matrix 



F = 



w° 


w° 


w° . 


w° 


w° 


w l 


w 2 . 


. W N ~ l 


w° 


w 2 


w 4 . 


jy2(JV-l) 




W N-l 


V^2(W-1) 





(13) 



where W k - Qxp(-2njk/N) and N = 16. The input data are the 
16 voltages measured at the lags, written as a column vector vo. 
The recovered spectrum is given by 



s = Fv . 



(14) 



Only the first eight of the elements in So are used because 
the remaining eight are related to the first eight elements (see 
Sec. 12.21) . This opens up a way for speeding up the DFT by 
adopting an 8 by 16 matrix. To keep things simple, we proceed 
with the full 16 by 16 matrix. 
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DFT usually assumes uniform sampling but Sec. 14.21 
showed that the sampling may be irregular by up to 10% be- 
cause of errors in the lag spacing. In the absence of noise, the 
true spectrum can be re covered from non-unifo rmly sampled 
data - see, for example Ba gchi & Mitra (1998). We have de- 
veloped a method that is particularly well suited for interfer- 
ometers and uses an astronomical source for calibration. It can 
be shown that the degree of correlated noise between the lags is 
minimised when Nyquist-sampled. In a correlator with lag er- 
rors, the noise becomes more correlated. The overall effect of 
non-uniformly sampled da ta is to degrade th e SNR when there 
is noise - see, for example IBracewelll ( 1999|) . The sampling er- 
rors are calibrated by defining a recovery matrix R that will 
return the true spectrum So, given the non-uniformly sampled 
lag data v; 



FRv = s . 



(15) 



R can be thought of as a transformation that interpolates the 
non-uniformly sampled data v to the uniformly sampled data 
vo that would have been measured had there been no lag errors. 
R also corrects for variations in the gain of individual detec- 
tors and dispersion in the microstrip delay lines. This is similar 
to t he scheme used for the WA. SP2 lag correlation spectrome- 
ter ( Harris & Zmuidzinasll2001 ). except that WASP2 was cali- 
brated using a CW signal to determine the recovery matrix R. 
Because it is a laborious task to calibrate each correlator indi- 
vidually, we use an unresolved astronomical source to calibrate 
all correlators simultaneously in situ. In an ideal correlator, the 
cross-correlation function is bound by a sine envelope. When 
Nyquist-sampled, the central peak of the sine envelope lies on 
a lag and the nulls of the envelopes coincide with the other 
lags (when r g = Ti). To form the 16 by 16 calibration matrix, 
the peak is aligned to the first lag. This will give a set of mea- 
surements like vi = (1,0.. .0) T , where T is the transpose. As 
the source is tracked, the peak drifts through in delay space and 
when it is aligned with the second lag, the second measurement 
V2 = (0, 1, . . . 0) T will be taken (Tig.[T4l). Repeating this pro- 
cess and gathering these 16 column vectors of measurements 
gives 



V = ( vi v 2 



Vl6 )• 



(16) 



In the ideal case, such as this, Vo is the identity matrix I. The 
matrix So is the corresponding collection of ideal spectra from 
the full calibration process. This means that So is equal to the 
Fourier transform matrix; 



FI = S = F. 



(17) 



In reality, there are gain variations, dispersion and non-uniform 
sampling. For small errors, V will be close to diagonal. We 
introduce the recovery matrix R to find the ideal spectrum So, 



FRV = S . 

Since S = F from (TJ). 
R = V" 1 . 



(18) 



(19) 



In the real system, the spectrum will not be rectangular. The 
spectrum for each detector may be recovered by taking the 



Fourier transform of the detectors' responses over the full cal- 
ibration run. The individual spectra of the detectors can be ac- 
counted for by gathering these measurements into the spectrum 
matrix S, then the recovery matrix is given by 



R = F SV~ 



(20) 



It can be seen from the DFT matrix (Eq. IT3l) that the first ele- 
ment of the spectrum vector So in Eqs. (ITU) and (TT3T) straddles 
the zero-frequency band. Because of the choice of lag spacings, 
this sub-band corresponds to 12 GHz IF. The information con- 
tent of this zero-frequency band is only half of what it should 
be because its two halves are redundant. At the other end of 
the spectrum, the highest frequency sub-band is split between 
positive and negative frequencies. This causes its centre fre- 
quency to be ill-defined. These undesirable symptoms can be 
avoided by shifting the sub-bands by half channel width us- 
ing the shift theorem for DFT: Multiplying the lag data by a 
set of complex factors shifts the sub-bands in frequency space. 
We define a diagonal matrix H whose diagonal elements are 
(6°, 6, 6 2 , • • • , e N ~ l ), where e = e~ nj/N . This channel shift matrix 
needs to operate before the DFT. The final recovered spectrum 
is given by 



s = FHRv. 



(21) 



Note that H should not be applied during the calibration step 
(Eq.[20|) because then H will be cancelled out. 



Lags 1 - 8 



Lags 9-16 




-0.2 -0.1 0.1 0.2 

HA /hrs 



-0.2 -0.1 0.1 0.2 
HA /hrs 



Fig. 14. Lag data during a calibration run on an astronomical 
source. To calibrate the correlator, a source is tracked with the 
path compensator held fixed. Unlike earlier plots which were 
plotted in instrument delay space (r g constant, Ti variable in 
Eqn[2]), these time- stream data plotted in hour angle which cor- 
responds to T{ constant, r g variable, and hence the fringe period 
is appropriate to the RF frequency (15 GHz) rather than the IF 
frequency (9 GHz). The matrix Vo (Eq.[l6| can be constructed 
from such a data set. 



10 



Christian M. Holler et al.: A 6-12 GHz Analogue Lag-Correlator for Radio Interferometry 



5.2. Simulations 

The calibration method was tested on simulated data with 
analytically-generated astronomical fringes assuming a flat 
spectrum. The spectrum was evaluated for an ideal correlator 
with no sampling errors and also with up to 10 % lag errors. 
Figures IT"5li and[T5b show that the calibrated spectra are sim- 
ilar to the true spectra, while the uncalibrated spectra deviate 
noticeably. The full spectra, including the negative frequencies, 
have been shown here. The two halves of the spectrum with 
channel shift (Fig.[T5b) are symmetric and redundant. 

The recovered spectra shown in Fig. [15] are not flat, de- 
spite the flat band-limited passband that went into modeling the 
cross-correlation function. Furthermore, the ripples across the 
spectrum change as the source drifts through the sky. This is il- 
lustrated in Fig. [16] Fig.[T6h shows a simulation of a processed 
data stream from a 5 m east- west baseline with a single point 
source at the centre of the field of view. After fringe-rotatior@ 
and standard phase calibration, we expect constant amplitudes 
and zero phase in all sub-bands. Instead we see periodic cycles 
in both amplitude and phase at the fringe rate for each sub- 
band. The magnitude of the cycles are worse for a correlator 
with lag errors. Calibration does reduce the magnitude of the 
cycles (Fig. [T6b ) but it remains noticeably higher than for a 
correlator with no lag errors (compare Fig. [T6b with Fig.[T6b). 

The cycles arise from aliasing. Because we can only sample 
the cross-correlation function over some finite range of delays, 
the recovered spectrum is a convolution of the true spectrum 
with the spectral response of the window function. But our 
signal is critically sampled, so the spectrum lies shoulder-to- 
shoulder with its aliased images (see Fig. [17]). The convolved 
wings of the images extend beyond the 6 GHz bandwidth and 
alias the signal's spectrum. This is why the edge sub-bands are 
most affected. Simulations show that narrowing the passband 
bandwidth by one sub-band reduces these alias cycles, as ex- 
pected. 

Applying a tapered window functions like Hamming win- 
dow before the DFT makes the alias cycles worse. The response 
function of the Hamming window has a wider lobe compared 
to a rectangular window (i.e. applying no windowing). So the 
recovered spectrum will have wider convolved wings and the 
alias cycles will be worse. However the Hamming window 
will suppress spectral leakage between channels because it has 
lower sidelobes (-43 dB) compared to a rectangular windows 
(sidelobes up to -13 dB). Applying the Hamming window re- 
duces alias cycles in the central sub-bands but the benefits are 
marginal (compare Figs.[T6h and[T8l). 

Window functions could be beneficial in correlators that 
oversample the cross-correlation function. For example, if the 
6-12 GHz signal is oversampled at 24 GHz , half the sub- 
bands between 0-6 GHz will have no signal. But these empty 



3 In interferometers with independently-mounted antennas, the 
phase centre drifts as the source is tracked. This changes the phase 
of each sub-band at a constant rate. The rate is determined by the ge- 
ometry of the telescope and the frequency of the sub-band. As we 
know this rate, the varying phase of the sub-bands can be made con- 
stant by a process called fringe-rotation. This is equivalent to stopping 
the fringes in delay space. 
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(b) With channel shift 

Fig. 15. Calibration can improve the recovered spectrum so that 
it is closer to what would have been measured by an ideal cor- 
relator with no lag errors. Here the lag errors are up to 10% 
of the lag spacings. (a) The top plot is the response of a point 
source in instrument delay space. The open circles are the volt- 
ages measured at each lag. The lower plots are the amplitudes 
and phases of the recovered spectra without and with calibra- 
tion. In (b), the channels were shifted. The effect of this is most 
dramatic in the phase. 
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(a) Rectangular window - with and without lag errors 
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(b) With and without calibration 



Fig. 16. Time-stream data of the amplitude and phase for each 
sub-band, (a) compares the time-stream data from an ideal cor- 
relator (dark curve) and a correlator with up to 10 percent lag 
errors (gray curve). It can be seen that lag errors makes alias cy- 
cles worse, (b) Calibrating the correlator with lag errors makes 
a small improvement to the magnitude of alias cycles. The gray 
curves are without calibration and the dark curves are with cal- 
ibration. The steps in the data are due to the path compensator 
changing the compensated path. 




v/GHz 



Fig. 17. A 6-12 GHz band-limited signal is sampled at = 
12 GHz. If we only sample over a finite range of the cross- 
correlation function, the estimated spectrum is no longer finite 
and band-limited. Although Nyquist's sampling theorem is ap- 
parently satisfied, the recovered spectrum suffers from aliasing. 
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Fig. 18. Time-stream data of the amplitude and phase for each 
sub-band like Fig.fT6h. This time, a Hamming window has been 
applied and the alias cycles are more severe, particularly in the 
edge sub-bands. The dark curves are for an ideal correlator and 
the gray curves are for a correlator with 10% lag errors. 



sub-bands will act as a buffer to protect the signal from be- 
ing aliased. A suitable window function with lower sidelobes 
will suppress aliasing further and reduce spectral leakage be- 
tween sub-bands. Applying a window function will necessar- 
ily broaden the effective spectral resolution. This will need to 
meet the upper bandwidth limit of the sub-bands (Eq. [3]). It is 
important to note that not all window functions are appropriate. 
The Hamming window will taper the sampled cross-correlation 
function. But it will also down- weight sources at the edge of the 
field of view because the peak of its correlation function will 
be offset from the central lag. An offse source's flux density 
will be systematically underestimated. This bias will vary over 
the course of an observation (see Fig. [19]), so it is impractical 
to correct for it. We have found that a raised trapezoidal win- 
dow function works well. It should be flat over the delay range 
of the field of view and tapered to non-zero levels at the edges 
(the edge lags contribute some useful information). Similarly, 
a raised version of the more curved Tukey window ca n work 
well. A review of window functions can be found in lHarrisI 
dl978h . 

The alias cycles are periodic, so the effects are reduced in 
the final data after smoothing. There will be some loss in the 
SNR but the introduction of systematic errors is a more serious 
concern. Residuals will be significant for long fringe periods 
comparable to the smoothing time. Parts of the data with long 
fringe periods will need to be edited out. The edge sub-bands 
suffer most from aliasing. In addition, the noise between sub- 
bands will be correlated. For these reasons, the edge sub-bands 
should either be rejected or down weighted. 

To summarise, our calibration scheme appears to recover 
the true spectrum. Simulations suggest that alias cycles will be 
present in the data as periodic cycles in the estimated spec- 
trum. These are caused by aliasing of the band-limited sig- 
nal. Applying a window function worsens the effect as does 
the presence of lag errors in the correlator. Calibration slightly 
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Fig. 19. The amplitude response for each frequency channel for 
a cross-correlation function centred away from the centre lag. 
This could arise from an off-centre source. The plot is the same 
as Fig. [3] except a Hamming window was used. The window 
will bias the measured flux density of off-centre sources. While 
alias cycles in the edge sub-bands (ch. 1 & 8) are increased, the 
alias cycles in the remaining sub-bands are slightly reduced. 
The tapered window also has the benefit of suppressing irreg- 
ular responses in Fig. [3] arising from sources at the edge of the 
field of view. 

reduced the magnitude of alias cycles. A software-based ap- 
proach that could overcome aliasing in a narrowband complex 
correlator will be discussed in a forthcoming paper. 

6. Manufacture and Commissioning 

The manufactured version of the correlators were packaged in a 
6U rack-mountable module that can be hot-swapped. The small 
signal from the detector needs to be buffered immediately, so 
the detector board was bonded to the amplifier board back-to- 
back. Bonding the PTFE board (with etching on the bottom 
side) to an FR4 board posed a significant manufacturing issue 
and could have been simplified had the slot-line phase- shifters 
been off-board. The whole assembly is housed in an aluminium 
box for shielding and the lid to the detector board is lined with 
microwave absorber to prevent standing waves. On the reverse 
side, a pair of op-amps amplify and remove the DC offset for 
each detector. The bandwidths of the amplifiers are matched to 
the video bandwidth of the detectors. 

6.1. Readout System 

A readout board slots next to the correlator module. It takes the 
16 signals from the correlator, demodulates and digitises them. 
Section 12.31 described how the IF signal at each telescope is 
phase- switch modulated by orthogonal functions / and g so 
that the cross-correlated signal can be extracted. The cross- 
correlated signal can be recovered by demodulating it with fg, 
which can be generated by taking the exclusive-OR of / and g. 
The readout board demodulates the signal by applying a gain 
of ±1 synchronously with the demodulation function. 

The number of bits needed by the analogue-to-digital con- 
verter (ADC) is determined by the required dynamic range 




(b) Amplifier side 

Fig. 20. The detector side of the correlator board (a) and the 
amplifiers on the other side (b). The correlator for each baseline 
is a self-contained unit. 



of the telescope. The dynamic range is given by the ratio of 
the minimum detectable correlated input signal and the corre- 
lated signal produced by the strongest source to be observed. 
At least three bits were assigned to recording the system noise 
and an additional 1 1 bits were needed for responding to bright 
sourcefl The readout system uses Texas Instrument's 20 bit 
ADC, DDC112U. The integration time of the readout system 
was chosen to sample the fastest fringes at least ten times in 
each fringe cycle. The effective sampling rate is 16 Hz but the 
data are stored at a lower rate where possible. 

The correlator system needs to be stable over a period of 
tens of seconds. The l/f noise of the correlator/readout system 
has a time constant of around 200 s with no IF input. When the 
IF inputs are connected, the 1 // noise is drowned by the tele- 
scope's system noise so the correlator system is dominated by 
the system noise. The performance was compared against the 
previous generation of readout system de veloped for the Ver y 
Small Array (a CMB interferometer, see I Watson et al J 12003 ). 
We were able to confirm that their performances were compa- 
rable, with the new system achieving slightly better SNR. 



4 We required sufficient dynamic range in the readout system to de- 
tect a 100 Jy source without saturating. We assume a telescope with a 
system temperature of 30 K and an effective area of 130 m 2 , together 
with an integration period of 1 /32 s. 
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6.2. Passband Response 

A problem with the passband of the correlator came to light 
after it had gone into production. The passband of the corre- 
lator was measured by feeding a swept CW signal into one of 
the ports. The CW was amplitude-modulated by a square wave 
and the correlator output was demodulated with a similar cir- 
cuit to the one used in the readout. Figure UJli shows a dip 
at 6 GHz and a fall off above 10.5 GHz. Compared with the 
original design (Fig. [9]), we had added a pin to take the signal 
to the amplifiers on the reverse side. The subsequent reactive 
load placed on the detector circuit led to the fall-off and coarse 
structures in the passbands (Fig.l2Th). 

We manufactured some individual detectors for use as to- 
tal power detectors to estimate the system temperature. These 
detectors used the same design as the correlators, except with a 
resistor in series with the diode to separate the transmission line 
from the video band circuitry. The passbands of this design are 
much flatter and do not roll-off below 12 GHz (Fig. EE})). The 
dip at 6 GHz is still present but at a smaller level. We believe 
that the 6 GHz dip results from the self-resonant frequency of 
the lumped components in the detector being comparable to 
our passband frequency. Smaller bulk components with higher 
self-resonant frequencies are now available and these may be a 
solution. 



7. Conclusion 

We have described the design and development of a broad-band 
(6-12 GHz) analogue lag correlator for radio astronomy appli- 
cations. Two design approaches were pursued: 1) the real cor- 
relator samples the cross-correlation function at Nyquist rate. 
2) the complex correlator samples both the real and imaginary 
cross-correlation function but at half-Nyquist rate. Both types 
were prototyped and tested. The 180° and 90° phase shifters 
used in the designs had errors less than ±6° and ±15° re- 
spectively over 5-13 GHz. There was very little difference be- 
tween the two designs in terms of performance. We adopted 
the real correlator for practical reasons. The passband of the 
manufactured version rolled off above 10.5 GHz due to a late- 
stage change to the design. We have shown that with a small 
modification, the full 6 GHz bandwidth can be achieved. We 
found unexpected errors in the lag spacings of up to 10%. We 
note that dispersion in the microstrip delay lines (±25° across 
the band) may cause problems, particularly for the outer lags. 
However, the effects are negligible in comparison to the lag er- 
rors. We have described a practical calibration scheme that uses 
an astronomical source to overcome delay errors and gain vari- 
ations in the lags, as well as dispersion. Simulations show that 
the edge channels will be aliased in a critically- sampled cor- 
relator like ours. We will explore this issue in a forthcoming 
paper. 
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Fig. 21. (a) The IF response of each lag over 5-13 GHz . The 
boundaries for each of the 0.75 GHz sub-band are shown by 
dotted lines. Note that unlike the passband for the earlier pro- 
totype in Fig. [T3k the passband drops off for the production 
version drops off above 10.5 GHz. (b) The IF response of 16 
modified detectors. 
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